On the possibility that STS "gap-maps" of cuprate single crystals are dominated by k 
space anisotropy and not by nano-scale inhomogeneity. 
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The results of a computer analysis of a simple 2D quantum mechanical tunnelling model are 
reported. These suggest that the spatial dependence of the superconducting energy gap observed 
by Scanning Tunnelling Spectroscopy (STS) studies of single crystals of the high T c superconductor 
Bi 2 Sr2CaCu2 08+cc is not necessarily caused by nanoscale inhomogeneity. Instead the spatial depen- 
dence of the STS data could arise from the momentum (k) dependence of the energy gap, which is 
a defining feature of a d-wave superconductor. It is possible that this viewpoint could be exploited 
to obtain k dependent information from STS studies. 
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I. INTRODUCTION 

In recent scanning tunnelling spectroscopy (STS) ex- 
periments, differential conductance (G = dl/dV) ver- 
sus voltage (V) curves are measured at many thousands 
of points on the surface of single crystals of the highly 
anisotropic cuprate superconductor E^S^CaC^Og+x 
or Bi:2212 1,2 . The form of G(V) at every point is roughly 
that expected from a tunnel junction between a normal 
metal and a <i-wave superconductor. However, the super- 
conducting gap parameter, A, determined from the peaks 
in G(V), varies strongly and continuously with position, 
from 20 to 70 meV with a typical length scale of 1 — 2 
nm. It is thought that G(V) is a direct measure of the 
local quasi-particle density of states (DOS) and therefore 
these "gap maps" are widely accepted as being associated 
with nanoscale inhomogeneity, implying that the super- 
conducting gap and all other superconducting properties 
are spatially inhomogeneous. It is difficult to reconcile 
this picture with the results of other experiments such 
as NMR 3 and heat capacity 4 , where such gross inhomo- 
geneity is not detected. 

The energy (E) dependent quasi-particle DOS, N(E) 
in a d-wave superconductor with a k-dependent order pa- 
rameter A(k) can be thought of as the sum of s-wave like 
contributions from different directions in k-space, i.e., 
N W = jfa I 4 S k)| [g2 _ A f k)2]1/2 > where dS is a small 
element of a constant energy surface in the normal state 
with electron velocity v{k). The measured quantity G(V) 



is proportional to 
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where P(k) is the probability of tunnelling into a state 
k, so G(V) only reflects the behavior of N(E) precisely 
when P(k) is constant. For a defect-free planar metal- 
insulator-superconductor tunnel junction, transverse mo- 
mentum (kr) will be conserved because of translational 
invar iance. Also, because of the exponential decay of 
the evanescent wave in the barrier, the tunnelling cur- 
rent will be dominated by states with smaller values of 
|kx|, that lie within a "tunnelling cone" whose angular 
width depends on the strength and width of the barrier—. 
So, for planar tunnelling into a superconductor with a 
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reasonably large Fermi surface, P(k) is not constant, be- 
cause final states with larger values of |kx| will not be 
accessible. But for planar tunnelling perpendicular to the 
conducting layers of a quasi-2D material such as Bi:2212, 
P(k) is constant if the quasi-2D Fermi surface is electron- 
like and exactly circular 6 , or if the "tunnelling cone" is 
wide enough, i.e. |ky| is large enough, so that P does 
not vary much over the Fermi surface 6 . At first sight, 
spatial confinement of the initial electron states in the 
sharp STS tip would be expected to give a larger spread 
in k^, so there would be an even greater tendency for 
P to be constant for c-axis tunnelling into Bi:2212 using 
STS. This seems to be the prevailing view in the field. 
While the model calculations reported here do show that 
for realistic STS tip sizes and tip materials, k^ is not con- 
served, i.e. there is no well-defined refracted wave in the 
quasi-2D metal, they also show that P is not constant. 
Instead it falls drastically as the in-plane wave number of 
Bi:2212, represented by the parameter \k y \ in the model, 
is increased. 

Recent STS studies give clear experimental evidence 
that bigger gaps are observed in areas of the crystal where 
there are more off-stoichiometric oxygen atoms in the 
BiO-SrO layer—. In a naive semi-classical picture these 
O 2- ions would scatter the tunnelling electrons, via the 
Coulomb interaction, giving them extra transverse mo- 
mentum. So there could be a tendency for a particular 
in-plane direction of k to have extra weight in G(V) at a 
certain separation between the STS tip and an O 2- ion 
and hence a particular value of A(k) could be favored 
there. This provides a mechanism for variations in the 
peaks in G(V) as the STS tip is moved across the surface 
of the crystal without invoking nanoscale inhomogeneity. 

The model calculations reported here are consistent 
with this semi-classical picture in that the strong k y 
dependence of P mentioned above becomes much less 
marked when there are O 2- ions near the STS tip. Al- 
though they do not give a completely adequate descrip- 
tion of the observed G(V) curves, they seem to capture 
the important physics and future microscopic calcula- 
tions of the G(V) curves should take the present results 
into account. 
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FIG. 1: (a) Spatial regions used for 2D tunnelling calcula- 
tions, I STS tip, II vacuum, III BiO-SrO layer, IV and IV 
oxygen ions, V Bi:2212. (b) Corresponding potential ener- 
gies versus distance (x) at y = 0. (c) Typical Fermi surface 
cross-section^, with solid arrows marking the minimum and 
maximum k vectors considered here and dashed arrows show- 
ing the angle a with the (— 7r,0) direction (see text). 



II. THE MODEL 

The spatial dimensions and energy barriers of the 2D 
quantum mechanical tunnelling model are defined in Fig- 
ures 1(a) and 1(b) respectively. Region I corresponds to 
the STS tip, which is the projection of a truncated cone 
with sides at 30 degrees to the surface normal, i.e. to the 
c-axis of the cuprate crystal. The Fermi energy of 5.5 eV 
in the tip corresponds to a free electron wave vector of 12 
nm which is close to the Fermi wave vector (k F Au ) of 
typical tip materials such as Au or a Pt alloy. The vac- 
uum gap in region II was generally taken to have a width 
La = 0.24 nm and a barrier height of 3 eV as deduced 
from STS work 8 . The boundary condition on the small 
line of length 2Lt represents an incoming electron wave 
of the form exp k J ^ u (x cos# + ys'mO)] travelling at an 
angle to the normal, i.e. on this line (where x=0) the 
real and imaginary parts of the wave function, i(jr and 
ipi were set equal to cos(kp u y smO) and sin(kp u y sin#) 
respectively. The STS tip width 2Lt was often fixed at 

0. 3 nm but calculations were made for a range of val- 
ues between 0.15 and 5 nm. For 2Lt = 5 nm (i.e. 10 
times the wavelength of the electrons in the Au tip) the 
result expected for a planar junction, namely conserva- 
tion of transverse momentum was obtained. However for 
2Lt < 0.6 nm there was no evidence for conservation, 

1. e. no sign of a wave propagating at the appropriate an- 
gle in region V. The angle 6 was varied between and 30 
degrees, but since the number of initial states in the tip 
increases as sin#, Ref. 5, it was decided that 25 degrees 
was appropriate for comparison with experiment. This 
choice does not affect our conclusions. 

Region III corresponds to the BiO-SrO layer which is 





0.7 nm 



FIG. 2: (Color on-line). Contour plots in the (x,y) plane of 
real part of wave function i/jr in region V for wavelengths of 
(a) 1.4 nm and (b) 0.7 nm, using the parameters given in Fig. 
1(b) and without any O 2- ions. Here x is measured from the 
start of region V. The linear color scale is also shown. 



taken to be 0.56 nm wide 9 . In most of the calculations a 
barrier height of 1 eV, consistent with band theory 7 , was 
used . This layer is now believed to be insulating in view 
of (a) the extremely large and non-metallic c-axis resis- 
tivity of Bi:2212 crystals^and (b) various experiments on 
intrinsic c-axis tunnelling in mesa structures, for example 
Ref. [Tl|. Region V, where the potential energy V(x,y) is 
constant, corresponds to the bulk of the cuprate super- 
conductor. Representing the electronic structure of the 
cuprates by a free electron model, which is effectively 
done here, is a drastic approximation. However includ- 
ing an extra barrier to represent the c-axis lattice poten- 
tial, i.e. the next BiO-SrO layer, does not alter the main 
conclusions because a defect-free planar barrier would re- 
flect semi-circular waves such as those shown in Figure 
2 without altering k y . This point was checked explic- 
itly by introducing an appropriate extra planar barrier 
in region V. This had no effect on the behavior shown 
later for P(k y ) in the absence of O 2- ions, but because 
of standing waves there was more scatter in the results. 

The present model is therefore a useful starting point 
for discussing the effect of oxygen ions in the barrier layer 
and for discussing whether the STS technique does indeed 
give a true k-space average of the electronic DOS. The 
value of the Fermi energy Ep in region V was varied 
over the range 0.76 to 3.04 eV because this corresponds 
to free electron wave numbers from 4.47 to 8.94 nm -1 
and wavelengths between 1.4 and 0.7 nm. These wave 
numbers are approximately equal to the minimum and 
maximum in-plane values of kp for electron states at the 
Fermi energy in the cuprates shown in Figure 1(c). In 
tunnelling problems one always considers the "physical 
electron" as emphasized by Anderson in connection with 
inter-layer tunnelling in the cuprates. 

The small circular regions IV and IV in Figure 1(a) 
correspond to spherically symmetric O 2- ions with a full 
outer shell situated 0.1 nm from the Cu02 layer. They 
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FIG. 3: Plots of the electron wave function (ipR, tyi and |^|) in a narrow strip of width x = 0.05 nm extending from y = -3.8 
to 3.8 nm, just inside region V, i.e. the first CuC>2 layer. Computed data for the two limiting wavelengths are shown in (a) and 
(b) in the absence of O 2- ions and in (c) and (d) with O 2- ions at +0.2 and -0.6 nm from the center of the STS tip. Parameter 
values are those in Figure 1 (b), with 2Lt = 0.3 nm, La=0.24 nm and (see text) = 25 degrees. 



give rise to a Coulomb field outside a typical atomic 
radius of 0.06 nm. In the metallic region V, this field 
will be strongly screened by the mobile electrons (on a 
length scale of order k^ 1 ). As often done in tunnelling 
problems^ this screening is represented by an equal and 
opposite image charge in region V giving a dipole field 
from each non-stoichiometric O 2- ion in the BiO-SrO 
layer. For simplicity, dielectric screening in the BiO-SrO 
layer is ignored since positive polarization charges at the 
II-III boundary can then be neglected. However in or- 
der to compare with experiments, where the mean spac- 
ing between O 2- ions is typically 0.8 nm, the Coulomb 
field was multiplied by a Gaussian attenuation factor, 
exp(— {y — yi) 2 / (J 2 ) where y\ is the y co-ordinate of the 
O 2- ion and a = 0.4 nm. This was done because the 
tunnelling electrons will be scattered by local increases 
in the Coulomb potential above a uniform background 
potential given by more distant O 2- ions. The Coulomb 
potential inside region IV was set equal to a constant, 
namely its value at the III-IV boundary. For simplicity 
it was also assumed that the dipole electric field does 
not extend into the vacuum region II or into the metallic 
region V. 



III. METHOD AND RESULTS 

A commercially available partial differential equa- 
tion (PDE) PC package 12 was used to solve the 2D 
Schrodinger equations for ifjR and ipi with the above po- 
tentials. An example program attributed to Backstromi 2 - 
and initial tests for the textbook problem of tunnelling 
through a ID potential barrier— were helpful in setting 
the boundary conditions (BCs), details of which are given 
in footnote ll4|. The BCs on the semicircular bounding 
surfaces shown in Figure 1(a) were chosen so that out- 
wardly propagating waves are favored and the absence 
of significant standing waves was verified by comparing 
plots of | -0| and along various lines in region V. Most 
of the calculations were done for the geometry shown in 
Figure 1(a) with L x = L y = 4 nm but checks were made 
for larger values of L X1 for a rectangular bounding re- 
gion and also for a cylindrically symmetric 3D case. The 
graphical output given by the PDE program was typi- 
cally in the form of contour plots of ipR such as those 
in Figures 2(a) and (b) or plots of i/jr and log |^| along 
certain lines (not shown). Arrays of 101 x 101 i/jr and 
ipi values over various user-defined regions of Figure 1(a) 
were also generated and these were processed further us- 
ing other commercially available PC packages. 

Figures 2(a) and (b) show examples of the contour 
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plots for the two extreme hp values mentioned above in 
the absence of any impurities in the barrier. Note that 
the shorter wavelength waves are more concentrated in 
the forward direction. As justified in detail later, this 
leads to the key result of the present work, namely P(k y ) 
is not constant for tunnelling between a sharp STS tip 
and Bi:2212 and therefore the G(V) curves will give a 
distorted image of the DOS. The barriers are relatively 
weak so that the width of the "source" along y just inside 
region V is considerably wider than that of the STS tip 
(0.3 nm). For the parameters given in Figure 1(b) and 
used in Figure 2, the regions of high intensity below the 
STS tip and 0.2 nm inside the conducting region V have 
FWHM in |^|, of Sy = 1.16 nm and 0.84 nm, for A = 
1.4 and 0.7 nm respectively. The ratios X/Sy differ by a 
factor 1.5 and this, together with the fact that X ~ 5y 
and 5x, must be be responsible for difference in spreading 
out of the waves shown Figures 2(a) and 2(b). Another 
important point is that, as mentioned in the previous 
section, there is little asymmetry between positive and 
negative y values, despite the fact that the wave coming 
in to the STS tip is at 25 degrees to the surface nor- 
mal, and therefore has significant transverse momentum 
k y — kp u sin(25) = 5.1 nm -1 along y. In contrast to the 
planar tunnelling case, no refracted wave that conserves 
k y is produced because the STS tip width, 2Lt = 0.3 nm 
is too small, as shown later in Figure 4(c). 



Computed data are shown in more detail in Figures 
3(a) to 3(d) as plots of i/jr, ipj and log 1^1 versus y at 
a mean distance along x of 0.25 nm into region V. The 
data shown are in fact the sum of 100 ip values between 
x = and 0.5 nm, now measured from the beginning 
of region V, for fixed values of y. A strip of this width 
corresponds to the first Cu02 bi-layer in Bi:2212, 9 but 
the same results are obtained for different strips e.g. be- 
tween x = and 0.3 nm. The plots are shown on 
a semi-logarithmic scale in order to see the behavior at 
larger values of y. There are no significant oscillations in 
|^| showing that the BCs do indeed produce propagating 
and not standing waves. The periods of the oscillations in 
ipn and ipi along y are very close to the two wavelengths 
used, 0.7 and 1.4 nm. Figures 3(a) and 3(b) show data 
with no O 2- ions in the BiO-SrO layer. Comparison of 
the two figures shows that the shorter wavelength has 
lower amplitude oscillations, by at least a factor 10, over 
most of the range of y confirming the differences shown 
in Figures 2(a) and 2(b). Figures 3(c) and 3(d) show 
the effect of having one O 2- ion at y = +0.2 nm and 
another at y = —0.6 nm, a spacing typically observed in 
the experiments 2 , the center of the STS tip being at y 
= 0. The key point is that in the presence of O 2- ions 
the amplitude is much less dependent on the wavelength. 
This suggests that STS does not give a true DOS aver- 
age for quasi-2D superconductors and that the presence 
of O 2- ions allows states with larger values of k y to be 
accessed. 



IV. COMPARISON WITH EXPERIMENT 

In order make an initial comparison with experiment 
the following procedure was used. The ip(y) data such 
as those shown in Figures 3(a) to 3(d) were multiplied 
by exp(— i2iry/X) or exp(z27n//A) and integrated over y 
to give the (complex) numbers a(k y ) and a(—k y ). These 
are the quantum mechanical overlap integrals between 
an initial electron state in the STS tip and plane waves 
of wavelength A propagating in the ±y directions in re- 
gion V. Hence the quantity \a(k y )\ 2 + \a(—k y )\ 2 is the 
required probability P(k y ) of electron tunnelling from 
the tip into a state in region V with \k y \ = 2tt/X. P(k y ) 
was calculated for a range of A values between 1.4 and 
0.7 nm. Since the Schrodinger equation is linear in ip 
and the solutions are obtained by matching ip and at 
the boundaries it is plausible, but not proved, that the 
same probability function P(k y ) would also apply to the 
physically realistic case of the Bi:2212 bands where, as 
shown in Figure 1(c), there is a range of |k| values with 
the same (Fermi) energy. Typical behavior of P(k y ) is 
shown in Figures 4(a) and 4(b) for various positions of 
the STS tip relative to one or two O 2- ions respectively. 
Generally P varies as Aexp(—B\k y \), where A and B 
are constants, but especially when the spacing of the two 
ions is ~ A, there is some extra curvature in plots of log P 
vs. k y which probably arises from two-beam interference 
effects. As shown by the curve ta in Figure 4(b), two 
nearby O 2- ions reduce P by 2-3 orders of magnitude, 
this is understandable because the height of the BiO-SrO 
barrier is increased locally from 1 eV to approximately 
5 eV. Figure 4(c) shows the behavior of P(k y ) for differ- 
ent STS tip widths (2Lt) ranging from 0.15 to 3 nm in 
the absence of O 2- ions. The exponential behavior men- 
tioned above persists up to 2Lt = 0.6 nm, after which 
there is a sharper fall associated with the tendency to- 
wards conservation of k y for a planar junction. 

These P(k y ) values were then used to calculate the cor- 
responding G(V) curves expected in the superconducting 
state. Referring back to Figure 1(c), the contribution to 
G(V) from a k- vector on the Fermi surface was weighted 
by P(|k|) where |k| is the distance from the Y point. Be- 
cause of the simplicity of the model, calculations were 
limited to the region between the two solid arrows in 
Figure 1(c). The behavior of P(k y ) for different O 2- en- 
vironments shown in Figures 4(a) and (b) means that 
small gap regions near the d-wave nodes are heavily fa- 
vored in the absence of O 2- ions but that scattering from 
O 2- ions allows higher gap regions to be accessed. Figure 
4(d) shows that this conclusion is not strongly dependent 
on parameter values such as the STS tip width 2Lt, the 
BiO-SrO barrier height or the spacing between the STS 
tip and the Bi:2212 surface (La)- Some typical calculated 
G(V) curves are shown as un-normalized plots in Figure 
5(a) and as normalized ones in Figure 5(b). In generating 
these curves the standard weak-coupling <i-wave form for 
the gap parameter, A(k) = Acos(2a) was used, where a 
is the angle between k and the (— 7r, 0) direction shown 



0.4 0.5 0.6 0.7 0.8 0.9 

k do'V) 



0.4 0.5 0.6 0.7 0.8 0.9 

kuo'V) 



^ 0.01 



1 1 1 1 1 | 


■■■(ill 


■ i i i i i 

1.5 tq 
3.0 tk " 






6.0 tn 


\ 


\ 




~— *— 12 to 


\\ 




18 tp 




^x0.01_ 


^^30tr 


\ 










(c) 

■ i i i i i i i i i 


x 0.001 


X 




-17 



P=Aexp(Bk ) 

/ (d) k 
i i i i 



0.4 0.5 0.6 0.7 0.8 0.9 

k (10 m ) 



1 2 3 4 5 
L(10" 10 m) or V(eV) 



FIG. 4: Plots of P(k y ), the probability weighting factor (see 
text) versus k y , (a) with no O 2- ions present (solid black 
circles) and with one ion present at different distances, in 
A° units (O.lnm) from the STS tip. (b) with pairs of ions 
at different distances, in A°, and (c) with no ions present 
but for various tip widths, 2L^p in A°. (d) Effect of various 
parameters, the BiO-SrO barrier height in eV, the vacuum 
gap in A°, and the STS tip half- width on the mean 
slope of log e (P) vs. k y in the absence of O 2- ions. 



by dashed arrows in Figure 1(c). Tunnelling data invari- 
ably show some broadening of the G(V) curves that is 
ascribed to a finite quasi-particle lifetime (h/T). In the 
Dynes formula^ used by Wei et al£ E is replaced by 
E — zT and N(E) is given by the real part of the usual 
expression. This formula should not be applied when 
r ~ A^and therefore, in the present calculations, mod- 
erate damping, namely T = 0.1A(k) was assumed for all 
k. In future it might be possible to derive experimental 
values of T(E,k) by appropriate fits to the STS G(V) 
curves. 



V. DISCUSSION 

The un- normalized G(V) curves shown in Figure 5(a) 
are similar to some of the published experimental data, 
e.g. the first paper by Pan et al. 1T where there is clearly a 
strong reduction in the magnitude of G(V) in the regions 
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FIG. 5: (a) Calculated G(V) curves for a <i-wave DOS in the 
same (arbitrary) units plotted versus eV/A for no O 2- ions, 
curve tk, and pairs of O 2- ions at various distances, |/ in A°, 
from the STS tip at y = 0. (b) Corresponding G(V) curves 
normalized to unity at eV = 3A and displaced vertically by 
multiples of 0.4. 



where the gap appears to be larger. It also noteworthy 
that curves labelled tk, tg, tf and sz seem to show very 
similar behavior at low voltages a feature that is also 
present in the experimental data, for example Figure 3 
of Ref. yj. The peaks in G(V) vary from V = 0.7A/e 
to A/e . However even though P(k y ) varies drastically 
with k y when there is little scattering by the O 2- ions, 
the corresponding G(V) curves do appear to be more 
rounded than those observed experimentally, especially 
recently 2 ^. 

Within the present approach even higher k- space se- 
lectivity may be needed to account for the experimental 
STS data. There are several ways in which this might 
occur. Firstly in the present model the band-structure of 
the BiO-SrO layer is ignored. In a microscopic treatment 
one would consider tunnelling to arise from virtual exci- 
tations into the BiO-SrO conduction band. This band is 
likely to have a narrow energy width and hence for a given 
impact parameter the transverse momentum gained by 
scattering from an O 2- ion may well be better defined 
than in the present model. Secondly the band structure 
of the Cu02 layer is not included. Since this is a tight 
binding band structure, Fermi surface states nearer the 
Brillouin zone boundary, i.e. those with small values of 
the angle a in Figure 1(c), will contain significant com- 
ponents of states with higher momenta, namely k + G 
where G is a reciprocal lattice vector. This will give even 
less weight to the states with higher values of |k|. 

Thirdly the standard "semiconductor model" is used 
to describe the tunnelling DOS of a superconductor. 
The theory of Blonder, Klapwijk and Tinkham (BTK) 19 , 
which makes use of the Bogoliubov equations rather than 
the Schrodinger equation, often provides good fits to 
G(V) curves of superconductors obtained with either 
metallic point contacts or small-area tunnel junctions, 
despite fact that it uses a ID model. BTK theory con- 
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tains some features that could be relevant here, for ex- 
ample the quasi-particle current transforms into a super- 
current over the coherence length which in the present 
case is highly k- dependent. A paper 20 applying BTK 
theory to a d-wave superconductor does conclude that 
G(V) rsj N(V), but because a ID model is used, k/r is 
conserved, in contrast to the present work. 

Fourthly, the presence of a pseudogap would also re- 
duce or eliminate contributions to G(V) from the regions 
of the Fermi surface with larger values of |k|. In a model 
used to account for heat capacity data 22 there is a tri- 
angular non-states-conserving pseudogap that ARPES 
experiments, e.g. Refs. [23| and 0, suggest has great- 
est effect in the anti-nodal directions. Introduction of 
a pseudogap would help in fitting the V-shaped, appar- 
ently non-states-conserving G(V) curves that are often 
observed experimentally by STS ^ 17 i 18 i 21 . Reports from 
other STS group a 21 i 25 i 26 suggesting that, especially for 
overdoped Bi:2212 crystals, the normalized G(V) curves 
can be uniform for linear scans over 10 — 20 nm do not 
necessarily contradict the picture presented here, since as 
shown in Figs. 4(a), 4(b) and 5(b) the normalized P(k y ) 
curves and hence the normalized G(V) curves, can be 
quite similar for a variety of O 2- spacings. 

Finally and perhaps most importantly, the k- 
selectivity could be further increased by the interference 
of scattered waves from different O 2- ions. There are 
some indications of this effect in the P(k y ) curve labelled 
tg in Figure 4(b) where the STS tip is placed symmet- 
rically between two ions at a distance ~ A from each. 
These effects could well be more marked in the realis- 
tic 3D case where "multiple beam" interference is more 
likely, for example the scattered waves from 3 or 4 suit- 
ably spaced O 2- ions could interfere constructively for 
certain directions of k to give 9 or 16 times higher in- 
tensity. Since the positions of the O 2- ions can be found 
using STS 2 it might even be possible to locate such struc- 
tures and hence obtain better k-resolution in the experi- 
ments. A faster PDE program would be needed to inves- 
tigate these 3D aspects using the present model. How- 
ever the main result reported here, namely the strong 
variation of P with k y in the absence of O 2- ions and 
the weaker dependence in their presence has been ver- 
ified for one special 3D case. This is the cylindrically 



symmetric situation where the O 2- ion is positioned im- 
mediately below the STS tip and the electron wave comes 
in at normal incidence. 



VI. SUMMARY AND CONCLUSIONS 

A new way of interpreting the "gap-maps" observed 
in STS studies of cuprate superconductors has been pro- 
posed that does not invoke nanoscale inhomogeneity. In 
principle, perhaps when extended to include real atomic 
orbitals and band structure, it can be used to obtain k- 
resolved information from STS data. 

Note added on 26th Jan. 2007. Since submitting this 
paper I became aware of analytical theory 27 of the scan- 
ning tunnelling microscope for a spherical STS tip. Equa- 
tions (4) and (9) of Ref. [271 lead to results that are con- 
sistent with the present work for a single barrier with- 
out any O 2- ions. Neglecting terms with reciprocal lat- 
tice vectors, ^ 0, gives a tunnelling probability P oc 
exp[-2(« 2 + |l?||| 2 ) 1/2 | -Tol], where 27 k = ^" 1 (2m0) 1 / 2 
is the minimum inverse decay length, k \ \ the in-plane 
electron wave-number, <j) the work function (i.e. the bar- 
rier height relative to Ep) and T^o the distance between 
the center of curvature of the STS tip and the sample 
surface. For a barrier height of 3 eV and \~r*o\ = 0.8 nm > 
fitting the above formula to P = Aexp (Bk y ) over the 
range of k y (= k \ \) used here gives B = -9.6 10 -10 m _1 . 
Bearing in mind that the present calculations were made 
for a flat STS tip and a 2D model, this agrees well with 
the data point in Figure 4(d), where B = -12 10 -10 m _1 
when the BiO-SrO and vacuum barrier heights are both 
3 eV, and the width of the combined barrier is 0.8 nm. 
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